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SUMMARY 

'tfk £ / 

A theoretical analysis is presented for obtaining hy use of 
Theodorsen 's propeller theory the load distribution along a propeller 
radius. to give the optimum propeller efficiency for any design condition. 
The efficiencies realized by designing for the optimum load distribution 
are given in graphs, and the optimum efficiency for any design condition 
may be read directly from the graph without any laborious calculations . 
Examples are included to illustrate the method of obtaining the optimum 
load distributions for both single -rotating and dual -rotating propellers. 


INTRODUCTION 


Recent contributions to the theory of propellers have been made by 
Theodorsen in a series of reports (references 1 to 4). In the first 
report of the series (reference l) a method based on electrical analogy 
was devised for obtaining the ideal circulation functions for single - 
rotating propellers. These circulation functions were shown to be in 
good agreement with the theoretical calculations made by Goldstein in 
reference 5 for two- and four -blade single -rotating propellers and with 
the extrapolations to other numbers of blades made by Lock and Yeatman 
in reference 6. The electrical -analogy method of measuring these 
functions was also applied to more difficult cases for which no theo- 
retical calculations had previously been madej in particular, to the case 
of dual -rotating propellers. 

Theodorsen in reference 1 introduced the concept of the mass 
coefficient k, which is an integrated value of the circulation functions. 
The mass coefficient represents the effective cross section of the column 
of the medium pushed by the propeller divided by the projected -propeller - 
wake area. 

This mass coefficient is made use of in the development of 
Theodorsen 's theory. In reference 4, expressions are given for computing 
the thrust, the energy loss, and the efficiency of any propeller with 
ideal circulation distribution based on the conditions in the final wake 
in terms of the mass coefficient. It is of interest to mention that the 
mass coefficient or mass of air operated on by the dual -rotating propeller 
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is much greater than that affected by the single -rotating propeller for 
the same set of operating conditions. This large difference in the mass 
coefficients for the two cases indicates that calculations for dual- 
rotating propellers based on the ideal circulation functions for single - 
rotating propellers are inadequate. 

Theodorsen 's theory, as previously mentioned, is based on the 
conditions in the final wake. The present analysis attempts to inter- 
relate the conditions in the final wake to the propeller and to give the 
information necessary to design a propeller for any desired operating 
condition- For single -rotating propellers, the method yields the same 
results as the conventional vortex theory with the Goldstein tip correc- 
tions applied. By the conventional vortex theory, however, it is 
necessary to determine the optimum blade -load distribution and then to 
make element strip -theory calculations in order to obtain the optimum 
efficiency for a given design condition. This procedure has been followed 
in reference 7 for a wide range of operating conditions. By Theodorsen 's 
theory the optimum efficiency t) can be obtained directly for any design 
condition from its relationship to the mass coefficient without laborious 
calculations. Thus, in the selection of a propeller for any design 
condition, a close estimate of the efficiency can be obtained before the 
design is made. 

The circulation functions and mass coefficients for the dual-rotating 
propeller were obtained in reference 1 for the ideal case and refer to 
conditions in the ultimate wake. Both propellers were assumed to operate 
in the same plane. Obviously, this condition is unattainable in the 
design of an actual propeller. The degree to which the ideal case can be 
realized in practice, or the applicability of the ideal functions to a 
given case, require further consideration and confirmation. 


SYMBOLS 


B 


number of propeller blades 


chord of propeller -blade element 


section drag coefficient 


section lift coefficient 


P. 


c 


ideal power coefficient (c s + e) 



total -power coefficient (P c + t r ) 


thrust 
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c Siji 

D 

d 

D o 

E 

Ej) 


e 


F 

K(x) 

l 

n 

P 

E 

r 

T 

t 

ta 

t r 

V 

v a 

V a 

Vi 


net thrust coefficient ( c s - t a ) 


diameter of propeller 
drag of propeller section 
diameter of wake helix surface 


ideal energy loss In wake ^ pFicw* 

energy loss due to blade drag 
induced energy loss coefficient 


(f v + l T )) 


E 


1 n V 3- 


pV J F 


projected area of helix (at infinity) 
circulation function 
lift of propeller section 

propeller rotational speed, revolutions per second 
input power to propeller 
tip radius 

radius to any blade element 
thrust of propeller 

power loss due to drag (nondimens ional) 
axial power loss due to drag 
rotational power loss due to drag 


forward axial velocity of propeller 
axial interference velocity (at propeller) 

average axial interference velocity (at propeller) 
resultant interference velocity (at propeller) 
rotational interference velocity (at propeller) 


k 
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V 
W S 
w 

V 
X 

a 

a i 

P 

X 


n ± 

p 

a 

ac l 

r 

<p 0 = tan 


average rotational interference velocity (’behind, each 
propeller) 

resultant velocity on the propeller at radius r 
local self -interference velocity 

rearward displacement velocity of helical vortex surface 

ratio of displacement velocity to forward velocity (w/V) 

radial location of blade element (r/R) 

angle of attack, degrees 

induced angle of attack, degrees 

blade angle, degrees 

advance ratio ( — - + — ) 

o J 

geometric advance ratio (v/jmD) 

/ F 

mass coefficient I 2 j K(x)x dx 

axial energy loss factor 

/ c s " t a\ 

propeller efficiency ( — j 

\ P c + t r / 

ideal propeller efficiency ^s/^c) 
mass density of air 

propeller element solidity (Bb/2nr) 
propeller element load coefficient 


circulation at radius x 




angle of resultant velocity W at plane of rotation 


-1 V/nD 


rtx 


co 


angular velocity 
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Subscripts: 

F 

R 

0-7 R 


front 

rear 

at 0.7 radius 


OPTIMUM PROPELLER DESIGN 
Single -Rotating Propellers 


Velocity diagram .- The velocity diagram for the single -rotating 
propeller is shown in figure 1. This figure is a reproduction of 
figure 13, reference 2, with some additional designations. The relation- 
ship between the axial interference velocity at the radius r , as given 
by the vortex theory, to the displacement velocity w of the vortex sheet 
is calculated in reference 2 and is shown in figure 1. The forward axial 
velocity of the propeller is V and the_tangential velocity with respect 
to the air at rest is cor. The vector bd is the resultant interference 
velocity V^ of the air with respect to the air at rest. Thus, the 

resultant velocity W of a point on the propeller at the radius r is 
given by the vector cd. The lift force l is perpendicular to this 
vector and the drag force d is exactly opposite in direction to W as 
indicated. From this figure a comparison of the method of analysis 
presented herein may be made with the conventional vortex -theory method. 

) It is required to find the point d in order to locate the end of the 
- velocity vector W and the angle cp that the vector W makes with 
- the direction of rotation. By the conventional vortex theory, the 
point d is located by starting with point b obtained from the V/nD 
of the undisturbed flow, proceeding in the V direction the distance V a , 


and then taking the perpendicular to this direction a distance V r . 

V + V 

(See reference 8 . ) The angle cp is given by tan cp = 


and W 


l + Va 

sin cp* 


cor - V r 

In the calculation of interference velocities V, 


and V r the local tip correction or Goldstein factor must be used to 
obtain the correct location of the point d. 


With the method developed in references 1 to 4, only the value 


of h v. 


which remains constant with radius, need he used. With this 


concept it is possible to use the integrated values of the mass coef- 
ficient as determined by the electrical analogy of reference 1 to obtain 
the detailed information needed at any radius. By this method the 


point d can be located by proceeding from point b 
the V direction to the point e and then down the direction of the 


a distance 1 w 
2 


m 
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velocity vector W a distance 
geometry of the figure as i w 


de , where 
sin 9 and 


de 


is obtained from the 


tan cp 


V + i w 
2 

a>r 


( 1 ) 


The resultant velocity is 

V + I w 

2 1 

W = — ; - — w sin 9 

sin qp 2 


= T” (V 

am cp \ 


i w cos 



( 2 ) 


The interference velocities may he obtained from the geometry of the 
figure by 


Vi = i- w cos 9 

1 2 

V a = Y ± cos <p = ^ w cos <p 


and 


V = v. 
r 


1 sin 9 = — w sin 9 cos 9 


Optimum blade -load distribution . - The design problem of an optimum 
propeller consists essentially in obtaining the value of the element 
load coefficient bc^ at each radius of the propeller blade. With the 
direction and magnitude of the relative velocity given at each station 
there remains only the choice of a section to give efficiently such a 
lift at the appropriate angle of attack. The value of c 1 should be 

at or near the ideal lift coefficient for the section in order to give 
minimum drag coefficient. 

The method developed in references 1 to 4 treats the velocity w 
as an independent parameter upon which all the other quantities depend. 
This reversal of procedure is convenient since all quantities are 
actually functions of w. The velocity w is related to the power 
coefficient P c of the propeller and also to the element load coef- 
ficient ac^ • The relation of w to ac\ is developed herein and the 
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relation of w to P c , which must he obtained in order to use it for 

design, is given subsequently in the section "PROCEDURE FOR PROPELLER 
DESIGN • " 

The required ideal circulation r(x) is given in reference 1 by 

r(x) = 2 » v)tf k(x) 

Bo 

= (UL^Ik(x) (3) 


In order to determine the element load coefficient bc^ the relation 

for the equality of the force on a vortex element and on an element of 
a lifting surface is given as 


pfW = — pW^c^b 


where b is the chord of the element. Hence, 


T = I Wc 7 b 
2 1 


W 


where W is given in equation (2), and thus 

1 1 A, 1 2 \ , 

r = o — : — T + rW COS cp ) c,b 

2 sin cp \ 2 / l 


(5) 


Using equations (3) and (5) for r gives at once the identity 


Do, = £L ±*h. K(x) 2 aln £ 

1 Bn 


1 2 
V + =• w cos cp 
2 


Introducing the nondimens ional velocity -w = the 
Bb V + i- w 

solidity a = - — •, and tan <p = (equation (l)) gives the non- 

2nr 2jtm 

dimensional relation 
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ac l 


( 


1 + w 




1 + 


i- V cos^p 
2 


2wK(x) 


. 2 

am 9 

cos 9 


( 6 ) 


The selection of a propeller for a given airplane installation may 
he based, on a method of evaluating a series of propellers for various 
operating conditions in order to determine the most suitable propeller. 

It is probable that several propellers, varying in diameter, blade 
number, propeller operational speed, and direction of rotation are 
equally as efficient for the design condition so that other considerations 
may enter into the propeller selection. However, the optimum efficiency 
for the propeller selected may be obtained from the charts, and therefore 
the load distribution along the radius that will give this optimum 
efficiency remains to be determined. 

The value of ac^ may be calculated for any radius from the 
relation 



1 + w 


+ 1 v 
2 


+ — w 

2 


2 

cos cp 


2wK(x) 


• 

sin cp 
cos 9 


where 


<P 


tan 


-1 


1 _V_ 

rt nD 





(7) 


Dual -Rotating Propellers 

In the design of dual -rotating propellers, it has been customary to 
select two propellers designed for single rotation and to use them as a 
dual -rotating propeller. • The fact that the circulation functions and the 
mass' coefficients obtained by the electrical -analogy method (reference l) 
are very much larger for the dual -rotating propeller than the sum of the 
values for the two single -rotating propellers indicates that the functions 
as used heretofore are not proper. The electrical -analogy method 
represents the case of, an idealized dual -rotating propeller in which the 
two components are in the same plane with the same load distribution on 
each component and with equal power absorption. Since actual propellers 
cannot conform to this ideal case, the applicability of the ideal 
functions requires further confirmation. Nevertheless, the optimum 
distribution for the dual-rotating propeller is essentially different 
from the single -rotating propeller, and in this analysis the loading 
functions and the mass coefficients as determined by the electrical- 
analogy method are assumed to apply to the optimum dual-rotating 
propeller. 
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Interference velocities for dual -rotating propellers . - The average 
axial interference velocity far "behind the propeller obtained from the 
momentum considerations is 


2V a = kw 


where k is the mass coefficient and w is the axial displacement 
velocity. This mean value is equally due to each of the two oppositely 
rotating propellers. The average axial interference velocity due to 
each is therefore exactly 



The average interference velocity at the propeller plane is one -half the 
value in the final wake and, therefore. 



where V a represents the average axial interference velocity at the 

propeller plane due to each component of the dual -rotating propeller. 
With the two propellers separated by a small axial distance, this 
velocity refers to a plane between the two propellers. The interference 
velocity at the front propeller is smaller and at the rear propeller is 
larger than at the plane between the propellers. In the following 
treatment, the propellers are considered to be very close together so 
that the axial interference velocity is the same on both propellers. 

In the final wake, the mean value of the rotational interference 
velocity for the ideal case is given by 

2V” r = 0 


For an infinite number of right and left blades equally loaded, rota- 
tional components would cancel exactly. However, the average rotational 
interference velocity immediately behind each propeller may be considered 
as 


V r = g- kw tan <p 
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In summary, the mean interference velocities acting on the front 
propeller from the rear propeller are: 


Axial: 



Eotational: 


= 0 


The mean interference velocities acting on the rear propeller from the 
front propeller are: 

Axial: — V „ = — kv 


Eotational: V r = i- kv tan 9 


It is useful to recognize that the mean self -interference of each 
propeller in its own plane is 


Axial: 


1 

4 


KW 


Eotational: 


i kw tan <p 
4 


Velocity diagram for the dual -rotating propellers . - The velocity 
diagram for the dual -rotating propellers is shown in figure 2. As in 
the case for the single -rotating propeller, the axial displacement 

velocity at the propeller is equal to i- v. In figure 2 the vector AB 
gives the mean axial interference velocity i kv of each propeller 

4 

acting on the other propeller. The vector BC gives the mean rotational 


interference velocity 


i- kw tan cp 


of the front propeller acting on the 


rear propeller. The total interference velocity acting on the front 
propeller from the rear propeller is therefore given hy AB, and the 
total interference velocity acting on the rear propeller from the front 
propeller is equal to the vector AC. The local self -interference 
velocity of the front propeller is given hy Wg , and the corresponding 


helix angle is given hy <£>• The local self -interference velocity of the 
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given "by cp^. The angle q>^ is slightly larger than the ideal helix 


smaller than cp. The design condition of most interest is the one for 
which T-p for each "blade of the front propeller is equal to for each 

"blade of the rear propeller. The number of "blades on the front and rear 
propeller are considered equal and the rotational speeds the same. This 1 
condition gives the self -interference velocity on the front propeller I 


equal to the self -interference velocity on the rear propeller and means 
that D and E must "be at the same horizontal level. 

As q>p and (p^ are needed in the design of the propeller, it is 

seen from figure 2 that the associated displacement velocity on the front 
and rear propellers has "been increased and decreased, respectively , "by 
the amount 




The displacement velocity is therefore 


Front: 



Rear: 



From figure 2, the velocity Wj> is shown to "be given "by the 
relationship 


Wn, = — + - kw sin cp, 

F sin <p 4 


o 


o 



o 



and the angle cp^, is given "by 
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tan 


< Pfp - 
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V + i-w^l + i-K tan 2 <p 


car 


V_ 1 _ 

nl) jcx 


i+i-vfi + i-K tan 2 cp 


(9) 


where <p Is given by the relationship 


V + i w 


tan qp = 


mr 


nD jcx 


Similarly, 


w r = - — + 7 kv sin <p_ + - kw tan <p_ cos ©^ 

sin © 4 0 2 0 0 


•V 3 

+ — kw sin © 

sin cp Q 4 0 


c 


,1 + t- kw sin 2 ® 
sin cp \ 4 o 


( 10 ) 


and 


tan <p = 

X\ 


V + 2 V ( 1 ' 2 K tan2<P ) 


car 


V_ 1_ 
nD jtx 


1 + ?r w (l - ir K tan 2 q> 

2 V 2 Y 


( 11 ) 


Optimum blade -load distribution * - The optimum distrihution of 
“blade loading is obtained from the determination of the element load 
coefficient bc^ at each radius from the fundamental relation 


i phc z W 2 = p rw 
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where f has been given in equation (3) by 


r 


(Y (• w)w 
Bn 


K(z) 


Eliminating r gives 


1 Bbc 7 W 

2 1 


(V + v)w 
n 


K(x) 


2 2nr ~ ° is tlle soli<il ' b y of each component of the dual -rotating 

propeller, if the number of blades in each component are assumed to be 
equal. Therefore, 


cc z W 


(1 + w)wV K(x) 
jrnDx v 


For the front propeller, this equation may be solved by use of 
equation (8) and 



X X 

nD itx 


(1 + w)v sin 9 


I p 

1 + — kw sin cp 

4 o 


iK(x) 


( 12 ) 


and for the rear propeller by use of equation (10) 

(”*) 


'E 


v , (1 + w)w sin <p 

- rr 2. K( x ) 

nD kx o _ 2 v 

1 + f - kw sin cp 
4 o 


(13) 


Use of Design Formulas 

In order to use the relation for crc^, note that it contains not 
only the independent variable w but also the function K(x) and the 

angle <p. The parameter K(x) should be expressed as a function of v + v 

nD * 

o 

which is based on the wake helix diameter. As was shown in reference 3 
however, D Q differs only slightly from the propeller diameter D and in 
the present design procedure . D is used instead of D Q . The function K(x) 


Ik 
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for single -rotating propellers is plotted against ^ + v in figures 3 

nD 

to 7* Similar plots for dual -rotating propellers were taken from 
reference 1 and are presented in figures 8 to 10- 


EQUATIONS FOE PEEFOBMANCE CALCULATIONS 
Single -Eotating Propellers 

In reference 4 the thrust has "been given "by 

T = pFKw 

and the ideal energy loss in the wake has teen given by 

E = pFKw 2 ^£w + 



With the introduction of the nondimens ional quantity w 
coefficient in nondimens ional form, is 


— , the thrust 
V 




and the induced loss coefficient is 


e - 


The power coefficient P c = c s 


E 



+ e is given ty 


( 14 ) 


( 15 ) 
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P, 


c 


2kv( 1 + w) ^ 


1 + — w 


K 


W ) 


( 16 ) 


The efficiency is given hy 



( 17 ) 


These formulas are all that are necessary for single -rotating 
propellers. The performance of the dual -rotating propeller is computed 
h y the same formulas . 


The coefficients c . e^ and P for the dual-rotating propellers 

o C 

are given in the same form as in equations ( 14 ), (15), and (l6) for single- 
rotating propellers. The only difference in the coefficients results from 
differences in the values of k, w, and which are substituted in 

the equations. 


Dual -Rotating Propellers 


The thrust of the front propeller is given by 



and with 



from equation (12) and Wjp from equation (8)' 



Similarly, for the rear propeller 


T J, = P S v 3 i? 



1 + — kw sin © 
v h 



dx (19) 
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Blade -Drag Losses 

The frictional loss or loss in efficiency due to the profile drag 
of the Blade is 

P 

= B g I bc d Vp dr 

Jo 


The drag force per unit length is h pW^bc,} where W, the resultant 

velocity of the "blade element, has been given in equation (2) for single - 
rotating propellers by 


V = 


sin <p 


0 


1 2 
+ - w cos <p 


) 


For the design condition, w is small, and because of the obvious 
uncertainties in the determination of the value of c d , it is not 

1 2 

necessary to retain the second term w cos cp. Introducing the 
solidity factor a = ^ - permits the drag loss to be given by 


Ep> = 


or in nondimens ional form 


2 3 

jtRpV ' 3 


ac d 


x dx 


|q sin <p 


E 


t = 


D 


\ pvW 


= 2 | -^1-x dx 

sin\> 


( 20 ) 


The component power losses are then, to the same degree of 
approximation in nondimens ional form. 


NACA EM No. L8F30 


17 


Eotational: t r 



ac a 3 , 

— x J dx 

sin 9 


( 21 ) 


Axial: 


t a = 2 


qc d 
sin cp 


x dx 


( 22 ) 


For the dual -rotating propeller operating at the design conditions, the 
terms containing -w are small, as is the case with the single -rotating 
propeller, and a close approximation to the drag loss is obtained if 
these terms are neglected. Furthermore, if it is assumed that the 
average of the resultant velocity W for the dual combination is equal 
to V for the single propeller, the equations (20) to (22) may be used 
for the dual -rotating propellers. Of course, for conditions other than 
the design condition, especially for very heavy loadings, exact drag- 
loss claculations require that the exact equations be used for either 
single -rotating or dual -rotating propellers. 

In su mm ary, the equations for obtaining the propeller performance 
are given by the quantities c s , e, and P c and the drag-loss factors 
are given by t r and t a . 

The net thrust power is 


T 


= c. 




V- 


(23) 


The power input is 




+ e + 


t r = P c 


+ t. 


(24) 


The efficiency is 


n 


P c + t r 



( 25 ) 


where from equation (1 6) 


P c = 2kw(1 + w) 
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The total power is also given by 


P = | P v\r 2 P Ct 


( 26 ) 


It should he remembered that the calculation is based on a given w* 
This procedure may seem unjustifiable since this parameter is not given 
by the specification but is the end result of a calculation based on the 
original data* The induced loss does not depend on the total -power 
coefficient P Cl p, but actually depends only on P c , and the quantity w 

cannot be obtained from the total -power coefficient. However, the 
/alue of P C| ^ in most cases exceeds P c by not more than 2 percent or 


p c = 0.98 Pc T 

Since P c in equation ( 16 ) is based on the w and the diameter 
of the final wake, and the value of P C| ^ in equation (24) is "based on 

the propeller diameter which is slightly larger than the diameter of 
the final wake, a very close approximation to w is usually given by 
equation (16). Therefore, 



2kw(1 + w) ( 1 + 


>( 



In some cases it may be necessary to calculate tr to obtain a more 
exact value of P c , especially if the blade profile drag is large* 


PROCEDURE FOR DESIGN OF PROPELLER 
Figures Used in Propeller Design 


The information necessary to design a propeller for any operating 
condition is given in the figures. Figures 3 to 7 give the circulation 
function K(x) interpolated for even fractions for two-, three-, four-, 
aix-, and eight-blade single -rotating propellers. The circulation 
function for the two -blade propeller was taken directly from reference 5 j 
lor the three -blade propeller, from reference 6 j and for the propellers 
having a greater number of blades was recalculated from the Goldstein 
tip correction factors as given in reference 7 . Figures 8, 9, and 10 
give K(x) for dual -rotating propellers with four, eight, and twelve 
blades, respectively. These values for the dual -rotating propellers 
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were taken from reference 1. Figure 11 gives the mass coefficient k 
for various numbers of blades for single -rotating propellers. Figure 12, 
which was taken from reference 1, gives k for dual -rotating propellers. 
The ideal efficiency t|^ is plotted against w for a range of 


in figure 15* The data for figures 13 and 14 were taken directly from 
reference 4 and the data for figure 15 was recalculated by the use of 
equation (l6) and figure 13* Figures 13 to 15 apply to either single - 
or dual -rotating propellers. The propeller efficiency may be calculated 
from either of these figures} however, in this report the efficiency has 
been determined from P c y4c as given in figure 15. 

Figures l6 and 17 give values of e, k, and c/k for two- and four- 
blade single -rotating propellers and figures l8 to 20 give values for 
four-, eight-, and twelve -blade dual -rotating propellers. The values 
of e for a propeller with a finite number of blades have not previously 
been published, but the values of e and c/k for an infinite number of 
blades are given in figure 4 of reference 4. The method for calculating 
e/* and e is given in the following section. 


In the selection of a propeller for a given airplane installation, 
the engine power, the forward speed, and the design altitude are usually 
specified. The selection consists of the determination of the number of 
blades, the propeller solidity, the propeller diameter, and the rotative 
speed. The ideal propeller efficiency, for any combinations of these 
variables can be readily obtained with the use of the charts . The 
procedure for a given blade number, propeller diameter, and rotative 
speed for either single or dual rotation is as follows: 

First, calculate the total -power coefficient. 



Propeller Selection 


P 


P 


2 4 


i pV3 2. D 2 


and then use this value for the ideal coefficient 





to find w 


• 1 
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It was shown in reference A that the dependence of the efficiency 
on g/ K in the efficiency formulas is very small and that it is suf f i - 
cient to know only the approximate value of e/ic* An examination of the 
formulas for c s and P c shows that their dependence on g/k is also 
small. It was further concluded in reference 4 that g/k is only 
slightly greater than k and that the practice of using c/k instead 
of k is considered satisfactory for design purposes. However, there 
appears to exist a simple relation between the axial -loss factor e and 
the total-loss factor k- This relation tales on the form of a 
differential equation 


_ i + 1 A. dK 

k “ 2 KdX 

where 


1 = 1 Y. t v 

* nl 0 


This relation has been checked and found to be exact for an infinite 
number of blades, and numerical checks for a two -blade propeller were 
in very close agreement. It is considered accurate for an empirical 
relation for design purposes for propellers of other numbers of blades. 

First obtain 


1 V_ 
S n nD 


as a first approximation to X. for use in the calculations. Then read 

off k and dn/dX from the appropriate charts of k against —(1 + w) 

nD 

for several values of w (figs. 11 and 12). Curves of e, k, and e /k 

are plotted against -L(i + w) in figures 1 6 to 20. Next plot a curve 

nD 

for the right side of the equation for P c against w- Where this 
curve intersects the horizontal line, P c = P c is the desired point. 

This value may be checked from the chart by inserting the values obtained 

from the plot in the equation. Thus are obtained k, w, —(1 + w) 

nD 

and g/k* From the chart of Pc/k (fig. 15 ), the optimum efficiency may 
be obtained. 
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The following examples illustrate the method of determining the 
optimum distribution of bc^, along the radius for both single -rotating 

and dual -rotating propellers that give the maximum possible efficiency 
(neglecting blade profile drag) that can be obtained with either 
propeller for one specified design condition. 


Illustrative Examples 


Single rotation . - Let the following data specify the propeller 
design conditions: 


Power, horsepower 2000 

Density, slugs per cubic foot O.OOIO65 

Velocity, miles per hour ^25 


The propeller selection has been made to the extent that the following 
data specify the propeller: 


Rotational speed, n, revolutions per second 

Propeller diameter, D, feet 

Number of blades, B 

V/n D 


The total P„ 


from the given conditions is 


. 23 

12 

k 

2.258 


'° T “ \ pV 3 rtR 2 


(2000) (550) 


|-(0.00l065)(623) 3 it(6) 2 


= 0.075 


The value of P c should be based on the wake diameter D Q instead 
of on the propeller diameter D and used to calculate w. Both P c and 
the contraction may be obtained by successive approximations but the two 
effects tend to cancel each other and generally P c ^ based on the propeller 

diameter is sufficiently accurate to use in the calculation of w. The 
relation between w and P c is given by equation (l6) as 


P c = 2Kw(l + w)^l + ^-w^ 
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where 


P 


c 



If values of v are selected to cover the range and the curve for the 
four -b lade propeller in figure 17 is used, the following table is 
obtained for the four-blade single -rotating propeller: 


w 

(assumed) 

K 

6 /k 

P c 

0 

0.245 

0.340 

0 

.1 

.215 

•313 

.0488 

.2 

.191 

.289 

.0970 


A plot of P c against w gives a value of w = O.I55 at P c = O.O75. 
Then, 


I-(l + w) = (2.258) (I.I55) = 2.61 

nD 

From figure 17, k is read at —(1 + w) =2.61, and the optimum 

nD 

propeller efficiency t]^ for a four-blade single -rotating propeller 
is read from figure 15* Thus 


k = 0.201 

P c 

— = 0.373 

K 


and 


Pi = O.929 

With w determined, ac^ for the single -rotating propeller may 
be found by a direct calculation from equation (6) 
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1 + v 


ac l = 


1 + | w 


, 1 - 2 
1 + 2 w cos <p 


27 K(x) 


sin^cp 
cos <p 


Values of the circulation function K(x) at each station are obtained 

from figure 4 at —(1 + 7) = 2.6l and the angle of the relative velocity 

nD 

at the propeller is given for each station by 


-> 1 - 

1 V P + 2 w 

tan <p = ± — 

jt nD z 


Performing these calculations for 7 = 0.155 gives the values of ac ^ 
and bc^ in the folloving table (the blade -width distribution, in feet, 
for a constant c^ of 0-5 is also given): 


X 

tan cp 

K(x) 

ac l 

bcj 

b 

(ft) 

A /0.7R 

0.1 

7.74 

0.033 

0-0842 

O.O79 

O.I58 

0.167 

.2 

3.870 

.078 

.0967 

.182 

• 364 

.386 

•3 

2.580 

.133 

.1054 

.298 

.596 

.631 

.4 

1-935 

.185 

.1044 

•393 

.786 

•833 

•5 

1.548 

.225 

.0952 

.449 

.898 

•952 

.6 

1.290 

.260 

.0855 

.483 

.966 

1.023 

•7 

1.106 

.271 

.0716 

.472 

.944 

1.000 

.8 

.968 

.257 

.0554 

.417 

.834 

.880 

•9 

.860 

.204 

• 0364 

•309 

• 6l8 

.655 

•95 

.815 

.146 

.0241 

.216 

.432 

.458 


Dual rotation . - The procedure is repeated for a 12 -foot -diameter 
four -blade dual -rotating propeller for the same design conditions as used 
for the single -rotating propeller. The folloving table is obtained for 
the four -blade dual-rotating propeller (values of k and e/ic were found 
from figure l8): 


V 

K 

«/* 

p c 

0 

0.472 

0.589 

0 

.1 

.432 

• 547 

.1002 

.2 

• 398 

.519 

.211 
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In this case a plot of P c against w gives a value of w = O.O75 
at P Q = 0 . 075 - Therefore , 


It is seen that the important parameter, the mass -flow coefficient, 
is 0 • kk 2 for the dual-rotating propeller and is only 0.201 for the 
single -rotating propeller. The effici'ency (without drag) is 9 6.1 percent 
for the dual -rotating propeller hut is only 92-9 percent for the single- 
rotating propeller. 

For the dual -rotating propeller the values of aci may he found 
for the front component from equation ( 12 ) } thus. 


— (1 + v) = 2.126 


K = 0-142 



K 


and 


q. = O.96I 



and for the rear propeller from equation (13) 



Equation (9) gives <Pp hy 
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and <p R is given in equation (ll) by 


tan 9 -r, = — — 
R nE jtx 


1 + — w 

( 1 - i K tan^q) ) 

2. 

V 2 V 


Values of the circulation function K(x) are obtained from figure 8 at 

the appropriate value of — — — = —(1 + w) . Performing these calculations 

nE nE 

for . V = 623 , n = 23 , E = 12, w = 0.075> and k = 0*442 gives the values 
of tan <p and crc^ in the following table: 


X 

K(x) 

tan qp F 

tan cp^ 

n 

Q 

O 

( tcj ) F 

(Nb 

0.1 

0-575 

IO.768 

4.145 

0.326 

0.321 

0.616 

0 .606 

•3 

.565 

2.608 

2.363 

.0995 

.0985 

.564 

•557 

.4 

.551 

I.916 

1.812 

.0692 

.0683 

.522 

.515 

•5 

•539 

1.518 

1.465 

.0501 

.0496 

.472 

.467 

.6 

•530 

I.258 

1.227 

.0370 

.0366 

.418 

.414 

•7 

.455 

I.O75 

I.O56 

.0268 

.0267 

•354 

•352 

.8 

•398 

•939 

.926 

.0191 

.0190 

.288 

.287 

•9 

•307 

•833 

.824 

.0122 

.0122 

.207 

.207 

•95 

•233 

— St 

•789 

.781 

.0085 

.0085 

.152 

.152 


A comparison of the optimum distribution of bc^ along the blade 

for the dual -rotating propeller from this table with the optimum distri- 
bution for the s ingle -rotating propeller as given in the preceding 
section shows that, if approximately constant c ^ is absorbed along the 

blade, wide differences in blade plan form will result for the two 
propellers designed for the same operating condition. For the operating 
conditions selected, the maximum bc^ for the single -rotating propeller 

occurs near the 0.6 radius and tapers rapidly towards the tip and the 
hub, being only slightly over l 6 percent of its maximum value at the 
0.1 radius. On the other hand, the minimum value of bc^ for optimum 

distribution for the dual -rotating propeller occurs at the propeller tip 
and progressively increases toward the inner radii. The value of bc^ 

at the 0.1 radius is four times its value at the O .95 radius. 

Since the design of the dual -rotating propeller calls for high 
loading over the inner sections, the efficiency of the dual-rotating 
propeller is less susceptible to compressibility losses which normally 
occur near the propeller tip for operation at high tip Mach numbers. 

The compressibility losses may be reduced by reducing the width of these 
sections or by reducing the operating lift coefficient. 
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Effect of "blade drag on efficiency . - The loss in efficiency due to 
the profile drag of the hlades can "be calculated from equations (20 ) 
to (22) if the "blade -width distribution and profile -drag coefficients 
at the operating c^ are known. Inasmuch as structural requirements may 

determine the shape of the "blade, especially over the inner radii, only 
one example is given. The equations, however, may "be applied to any plan 
form. The example selected is for the four -"blade single -rotating 
propeller on which the induced efficiency has "been previously calculated. 
The shank sections of the propeller "blade were assumed to "be round, 
similar to the Hamilton Standard Propeller No. 3155-6 and the "blade plan 
form from x = 0.3 to x = 1.0 was made optimum for a cj of 0.5. The 
profile -drag coefficients for the several radii are the same as given in 
reference 7 for the Hamilton Standard Propeller No. 3155-6 which has 
Clark Y sections and are given in the following table. It is assumed that 
a spinner covers the inner 0.2 of the radius. The distribution of oc l 

with x and of sin qp with x have been included in the table: 


X 

a ci 

a 

cd 

sin tp 

ac. 

& x - 

ac d x 3 

sin qp 

sin <p 

0.2 

0.0967 

0.1934 

0.400 

O.968 

0.01600 

0.00064 

•3 

.1054 

.2108 

.100 

.932 

.OO697 

.00061 

.4 

.1044 

.2088 

.020 

.889 

.00188 

.00030 

•5 

.0952 

.1904 

.010 

.840 

.00113 

.00028 

.6 

.0855 

.1710 

.008 

.790 

.00104 

.00037 

•7 

.0716 

.1432 

.007 

.742 

.00095 

.00046 

.8 

.0554 

.1108 

.006 

.696 

.OOO77 

.00049 

•9 

.0364 

.0728 

.006 

.652 

.00060 

•00049 


Performing the integrations and substituting in the formulas gives 
for rotational -drag -loss coefficient 


t r 


' 1.0 



x3 dx 

sin qp 


= ■ 2 . (0.000348) = 0.0014 
O.516 


and for the axial-drag -loss coefficient 




\ I MV »"* 
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' 1.0 


*a = 2 


' 0.2 


<* c d 
sin p 


x dx 


= 2(0.00213) = 0.0043 


The induced thrust coefficient has "been given by equation (14) as 


c s = 2kw 


1 + w( i + £ 

2 K 


= 2(0.201) (0.155) 


1 + 0.155 I o 7 0-2 9 


(I 


= 0.0700 


and the induced power coefficient hy equation (l6) as 


P c = 2 kw(1 + w) ( 1 + ^- w 


= 2(0-201) (O.I55) (1.155) (1-045) = 0.0754 


The induced efficiency is 


^ = pT 


0.0700 

0.0754 


= 0.929 


With drag included, the total thrust is given hy 


c s T “ c s t a 


= 0.0700 - 0-0043 = 0.0657 
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and 


p c T - p c + 


0.0754 + 0.0014 = 0.0768 


The efficiency is 



Thus it is seen 
preceding table 
to 85.5 percent 


that the "blade drag of the magnitude given in the 
reduces the propeller efficiency from 92*9 percent 
for the propeller operating conditions given- 


CONCLUDING EEMARKS 


A comparison of Theodorsen 's propeller theory with the conventional 
vortex theory shows that the optimum load distribution along the blade 
for single -rotating propellers obtained by the two theories is essentially 
identical and as a Result the optimum efficiencies- are the same for a 
given operating condition. Theodorsen's theory has the advantage, - however, 
that the optimum efficiency for any design condition can be obtained 
quickly and accurately by the use of the mass coefficient • k without any 
laborious calculations and before the final design is made. 

The distribution of the circulation function K(x) for the idealized 
dual -rotating propeller is radically different from the existing values 
for the single -rotating propeller that have been previously used for the 
dual -rotating propeller. Also, the mass coefficient k for the dual- 
rotating propeller is larger than the sum of the values for two single - 
rotating propellers. These quantities, which are not available from 
mathematical computations but are obtained from the electrical -analogy 
method of Theodorsen, are used herein for obtaining the optimum load 
distribution along the blade for the dual -rotating propeller. 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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Figure 1 .- Velocity diagram for single -rotating propeller. 
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Figure 2.- Velocity diagram for dual-rotating propeller. 
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nE ' nB (l + w) 


K(x) 


Figure 3.- Circulation function K(x) for two-blade single -rotating 
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Figure 4.- Circulation function K(x) for three-blade single -rotating 

rvronpllpr (rpfprpnop 
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Figure 5.- Circulation function K(x) for four -blade single- 
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Figure 6, 


Circulation function K(x) six-blade single- 
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1 


K(x) 



Figure 7.- Circulation function K(x) for eight -blade single - 

rotating propeller. 
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Figure 8.- Circulation function . K(x) for four-blade dual -rotating 
propeller (reference 1). 
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K(x) 



Figure 9.- Circulation function K(x) for eight-blade dual-rotating 
propeller (reference 1) . 
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_v_ 

nD 


(1 + w) 


K(x) .5 


Figure 10.- Circulation function K(x) for twelve-blade dual -rotating 
propeller (reference 1) . 











_ V 4- W 

Figure 11-' Mass coefficient K against for various numbers of blades 

nD 

for single-rotating propellers. 
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Figure 12,- Measured values of mass coefficient k for dual-rotating 
propellers with various numbers of blades. Two -blade single - 
rotating propeller included for comparison (reference 1) . 
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Figure 13.- Propeller efficiency against w (reference 4). 
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Figure 14.- Propeller efficiency against c s > (reference 4). 




























44 


NACA EM No. L8F3O 



Figure 15.- Propeller efficiency against P c /k. 
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Figure 16.- Values of € and k for two-blade single- 
rotating propeller. 





k6 


NACA EM Wo. L 8 F 3 O 



Figure 17.- Values of ? and k for four -blade single - 
rotating propeller. 




NACA EM Wo. L 8 F 3 O 


k7 



Y_J_w = JL (i + w) 
nD nD 


Figure 18.- Values of € and k for four-blade dual- 
rotating propeller. 
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Figure 19.- Values of e and k for eight-blade dual- 
rotating propeller. 
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Figure 20.- Values of € and * for twelve-blade dual- 
rotating propeller. 




